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Abstract 

Motivated by some new experimental data, we carry out a phenomenological analysis 
of D — > Kir decays including both Cabibbo favored and doubly Cabibbo suppressed 
modes. Two asymmetries, R(D°) and R(D + ), which are generated through interference 
between Cabbibo favored and doubly Cabibbo suppressed D — > Kir transitions, are 
predicted. The relative strong phase, 5k-k, between D° — > K~ir + and D° — ► K + ir~ 
decays, is estimated. The theoretical results agree well with the current measurements. 
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1 Introduction 



Non-leptonic D — > Xtt decays and their strong phases have been of great interest as they 
are essentially related to the studies of CP violation, D° — D° mixing, and SU(3) symmetry 
breaking effects in charm physics [H [21 El HI E] • These decay modes contain both Cabibbo fa- 
vored (CF) and doubly Cabibbo suppressed (DCS) transitions, and the effective Hamiltonian 
relevant for them is given by 




+ V us V* d [Ci(diCi) V ^A(UjSj)v-A + C 2 {diCj) V - A (UjSi) V -A]} 

+H.c, (1) 

where V — A denotes 7^(1 — 75). The first line in eq. (jTJ governs CF decays and the second 
line DCS decays. 

Theoretically, factorization hypothesis has been widely utilized in the hadronic D decays. 
Many studies are based on the naive factorization approach, which simply replaces the 
matrix elements of a four-fermion operator in a heavy-quark decay by the product of the 
matrix elements of two currents. This approach has long been used in phenomenological 
applications, although there is an obvious shortcoming that it cannot lead to the scale and 
scheme independence for the final physical amplitude. Several years ago, the authors of Ref. 
[6] have formed an interesting QCD factorization formula for the two-body exclusive non- 
leptonic B decays, in which the scale and scheme dependence of the hadronic matrix elements 
is recovered, and the naive factorization can be obtained as the lowest order approximation. 
The radiative corrections in the strong coupling constant a s , which are dominated by hard 
gluon exchange, can be calculated systematically using the perturbative QCD in the heavy 
quark limit. This means the strong final-state interaction phases, which arise from the hard- 
scattering kernel, are calculable from first principles. Analogously, in the heavy charm quark 
limit, a similar factorization formula for the matrix elements of the operators Q^s in the 
effective weak Hamiltonian ([1]) can be written as [6] 

(PiP2\Qi\D) = £ F^{ml) jf 1 duT^u)<p P2 {u) + (P 1 <- P 2 ) 

+ /' d£dudvTl\t,u,v)MO<i>pM<i>pM, (2) 

Jo 

where F* > ~* Pia (m\ 2 ) denotes a D — > P\ 2 form factor, 4>x(u) is the light-cone distribution 
amplitude of meson X. T-j(^,u,v) and Tf 7 (£, w, t>) are hard-scattering functions, which are 
perturbatively calculable. Then theoretical results for D decays can be obtained straight- 
forwardly. Taking the CF decays D° — > K~tt + , D° — > K°ir°, and D + — > K°tt + as examples, 
for the leading power contribution, we get 

B(D° -> K--K+) = 3.97%, B(D° K°n°) = 0.08%, B{D + -> R\ + ) = 7.66% (3) 

at the scale // = 1.5 GeV. (Here we have parameterized Jq 1 dt; <£d(£)/£ = m o/^D and set 
\d = 0.3 GeV in the numerical calculations.) The corresponding experimental data from [7] 
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are 

B(D° -> K--K+) = (3.80 ± 0.07)%, B(D° -> XV) = (2.28 ± 0.24)%, 



£(£>+ -> iTV + ) = (2.94 ± 0.12) 



/o. 



It is seen that, although the predicted branching ratio for the color- allowed decay D° — > 
i^ _ 7r + is in qualitative agreement with the data, the prediction for the color-suppressed 
decay D° — > ^°7r° is too small, and for the charged mode, the theoretical B(D + — > fT°7r + ) 
is too large. Similar conclusion will be reached when applying the formula ([2]) to singly 
Cabibbo-suppressed (SCS) and DCS decays. This seems to indicate that the charm quark 
mass is not heavy enough to apply the QCD factorization approach [6J or pQCD [5] in D 
decays very reliably. Therefore one generally appeals to the phenomenological analysis of 
these processes. 

Experimentally, many new results in D decays are expected soon from the dedicated 
experiments conducted at CLEO, E791, FOCUS, SELEX, and the B factories BaBar and 
Belle. In particular, as pointed out in [91 [10], there are interesting asymmetries due to 
interference between CF and DCS D — > Kit transitions, defined as 

mm = B{D - Ks7T) ~ B{D - Kl7i) (a) 

K ] ~ B(D -> K s ir) + B(D -> K L ir) ' { } 

which have been observed by CLEO Collaboration [TTJ very recently, 

R(D°) = 0.122 ± 0.024 ± 0.030, R(D + ) = 0.030 ± 0.023 ± 0.025. (5) 

Also a preliminary result on the relative strong phase between D° — > K~n + and D° — > 
K + 7i~, which is due to SU(3) symmetry breaking and important in the search for D°— D° 
mixing [HEI2], has been reported by CLEO Collaboration as 

cos 6 K7V = 1.09 ±0.66 PI, (6) 

although with very large uncertainty. 

Motivated by the new measurements mentioned above, we would like to perform a phe- 
nomenological analysis of both CF and DCS D — ► Kir decays. As will be shown below, the 
present data cannot allow us to determine all of the phenomenological parameters appear- 
ing in decay amplitudes. Implementing the SU(3) symmetry may constrain the amplitudes, 
thus largely reduce the number of independent parameters. However, it is known that this 
symmetry is not well respected in nature, even badly broken in some cases. Therefore, as 
a conservative way to constrain these amplitudes, in this paper we assume that SU(3) sym- 
metry in D — ► Kit decays is moderately broken, namely, symmetry breaking effects in decay 
amplitudes are dominated by decay constants fp and D — ► P (P = n, K) form factors, and 
other SU(3) symmetry breaking sources can be neglected. This is not a general feature of 
SU(3) symmetry breaking in charmed decays. 



2 D — >• Ktt decay amplitudes 

We begin by considering the D — > Kix decay amplitudes in terms of the quark-diagram 
topologies T (color-allowed), C (color-suppressed), £ (PU-exchange), and A (PU-annihilation) 
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[H], which are given by 



A(D° - 




= i|v^ + ^ 


(7) 


V2A(D° 


- /?°7r° 


)=i^WE(C-£), 


(8) 


A(D + - 


-> ^°7T + ) 


= i^^(T + C), 


(9) 


A(D° -> 




= <^VL.VS(T' + 0, 


(10) 


V2A(D° - 






(11) 


A(D + - 


► K°7T + ) 


= <^K.VS(C / + ^ / ), 


(12) 


V2A(D + - 






(13) 



and two isospin relations 

A(D° -> ir-vr + ) + v^A(D° -> K\° = A(D + -> i?°vr + ), (14) 

A(D° -> iT+TT-) + V^A(D° -> tfV) = A(D+ -> ^°7T + ) + V / 2A(D+ -> K + 7T°) (15) 

are satisfied explicitly. For our notations, we have extracted the CKM matrix elements 
and factor Gf/V% from the quark-diagram amplitudes, and the prime is added to DCS 
amplitudes. 

The present experimental status of D — > Kir decays is not very satisfying. Branching 
ratios of three CF modes and B(D° — ► K + tt~) have been reported by Particle data group 
[7j, however, there are no measurements for B(D° — ► fT 7r°) and B(D + — ► i^°7r + ) yet. Only 
an upper bound on B(D + — > i^ + 7r°) < 4.2 x 10~ 4 (CL=90%) is shown in [7]; while, very 
recently, BaBar Collaboration and CLEO Collaboration have given 

B(D + -»■ K+tt ) = (2.52 ± 0.47 ± 0.25 ± 0.08) x 10" 4 [E], 

(16) 

B(D + -> if+vr ) = (2.25 ± 0.36 ± 0.15 ± 0.07) x 10~ 4 pi], 

respectively. In general the quark-diagram amplitudes in (J7j) - ( FTBl could have non-trivial 
strong phases. Therefore, only using the available experimental data, it is impossible to 
determine these amplitudes without any theoretical assumptions. 

On the other hand, with the help of the factorization hypothesis, T, T 7 , C, and C can 
be expressed as 



T = 


U( m D - 


~m\) 


Ft 


*K 




af. 


C = 


fi<(m 2 D 


-ml) 


Ft 






af. 


T = 


fi<(m 2 D 


-ml) 


F D - 






af. 


C' = 


fK(m 2 D 


-ml 


)F D - 


->7T 


(m^ 
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(a) 



(b) 



(c) 



(d) 



Figure 1: Weak annihilation diagrams (W-exchange or W-annihilation) via gluon emission. 
The solid square denotes the weak vertex. 



where af ff 's are regarded as the effective Wilson coefficients fixed from the data (in the 
naive factorization, = Ci,2 + (^2,1/ H c ). Generally in the factorization approach af ff 's in 
DCS amplitudes could be different from the ones in CF amplitudes due to SU(3) symmetry- 
breaking effects. Here we do not distinguish them because, as stated above, we assume that 
the SU(3) symmetry breaking effects in D — > Ktx modes have been mostly captured by the 
decays constants (/„. and fx) and form factors (F D ~* n,K (q 2 )) in the amplitudes. 

A similar analysis of PF-exchange and ^-annihilation amplitudes leads to 



£ 
£' 



fD(m 2 K 

2 



m n a. 



efi 



f D {mt-m 2 K )F»- 
A' = f D (ml-m 2 K )F°- 



*K 



D) u '2 > 
eflf 



,cff 



{m D ) ai 



This will give 8 = —£'. However, the annihilation form factor F®~* (m 2 D ) is expected to 
be strongly suppressed due to the large q 2 = m 2 D [17], therefore contributions from eq. (jTSl 
are believed to be negligible. For B mesons, the weak annihilation amplitudes (H^-exchange 
or H^-annihilation) induced by the topologies of gluon emission arising from the quarks of 
the weak vertex, as shown in Fig. 1, have been analyzed in [18] . which is thought to be 
numerically important in B decays. The similar study on D mesons have been done in 
[T9| [20] . and it has been shown that these contributions could also play important roles in 
D decays. The 0(a s ) contribution can be read directly from Refs. 



£ 



Cf 



Hi 



c 



7T 



and 



A! 



H 



c 



18 \X A 



£' = £, 




(19) 



(20) 



where Xa = Jo dy/y has been used to parameterize the logarithmically divergent integrals 
due to the end-point singularity, and C\, Ci are the Wilson coefficients in ((T]). Note that the 
asymptotic form of the light-cone distribution amplitudes for light mesons have been used 
in the derivation of eq. (|T9l) . This is consistent with the assumption on SU(3) symmetry 
breaking used in deriving eq. (TT7|) . 

Meanwhile, comparing eq. fl20|) with eq. f|T9|) . one can get an additional constraint 

C 2 



A' 



Ci 



£'. 



(21) 
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A^ = 

MS 


215 MeV 


A^ = 

MS 


325 MeV 


A^ = 

MS 


435 MeV 


Ii [GeV] 


NDR 


HV 


NDR 


HV 


NDR 


HV 


1.0 


-0.339 


-0.390 


-0.400 


-0.464 


-0.464 


-0.548 


1.5 


-0.277 


-0.319 


-0.318 


-0.369 


-0.357 


-0.419 


2.0 


-0.240 


-0.278 


-0.271 


-0.316 


-0.301 


-0.353 



Table 1: The scale and scheme dependence of C 2 / C\ at the next-to- leading order. The values 
of C\ and C 2 are taken directly from (21]. 



We would like to give some remarks here. 

• By combining eq. f[2Tj) with eqs. f[T9l and f[TT|) . we will reduce independent complex 
phenomenological parameters appearing in decay amplitudes as af 1 , a| ff , and £. Since 
only five branching ratios of D — > Ktc decays are measured up to now, this means that 
including the additional constraint (I2TI) is important to enable us to determine the 
D — > K7r amplitudes from the present data. From [21 J, C\ and C 2 have opposite signs, 
therefore weak-exchange £ and weak-annihilation A amplitudes have opposite signs, 
consistent with the observations in [22]. Also |C 2 | < \Ci\ [2TJ, we have \£'\ > \A'\; 
while the contrary conclusion will be obtained if we use eq. ([TBI . 

• Strictly speaking, we have to admit that eq. (12Tj) is not very physical since C\ and C 2 
are both scale and scheme dependent [2TJ. The scale and scheme dependence of G 2 jC\ 
has been shown in Table 1, from which it is found that C2/C1 is about —0.5 ~ —0.3 
for around 1.0 ~ 1.5 GeV (Note that the scale in this range is relevant for D decays). 
Therefore, we will treat in the following numerical calculations C2/C1 as a negative 
parameter instead of a ratio of two Wilson coefficients. 

• The weak-annihilation contribution is power suppressed in the heavy quark limit. The 
divergent integral appearing in ffl9l and (TSUI) signals that factorization breaks down 
actually. In the present analysis we therefore use relations £ = £' and (12T1) for these 
three weak-annihilation amplitudes instead of their explicit expressions shown above. 
Although they are not model independent relations, one will find that phenomenolog- 
ically they work very well in D — > decays. 



3 Phenomenological analysis 

From now on we study some possible phenomenological applications based on the above 
theoretical assumptions. 
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First, the use of eqs. (fT7|) and (|T9j) gives C = C and £ = £', hence we have 

A(D° -> #V) = ^§A(D° -> #V) = - tan 2 9 C A(D° -> W), (22) 

which implies that the relative strong phase between these two amplitudes vanishes. Here 
6c is the Cabibbo angle. Consequently, one gets 

Using tan#c — 0.23, R(D°) ~ 0.106, which is in agreement with the measurement in eq. (J5]). 
The same result has been obtained in Refs. [9j [23]. However, one cannot expect the similar 
result as (|23|) for R(D + ) since there is no similar relation as (122!) between A(D + — > A'°7r + ) 
and — > fT°7r + ), even in the SU(3) symmetry limit. We will discuss this issue later. 

Second, using the constraint (|2T|) together with (|20|) and fflTj) . one can obtain an in- 
teresting relation among the amplitudes of D° —>■ K^tx^ and D + — > K + ir° decays, which 
is 

tan 2 6 C A(D° -> £T7r + ) + /^(D -> K+tT) = v / 2C^(^ + -> (24) 

where 

K ~ l + Ca/d ' ^ l + C 2 /C x 

with 

= K> - m^i^m 2 ) 

Thus the relative strong phase 5^ between — > i^ + 7r~) and A(D° — > i^~7r + ) is given 
by 

tan 4 6 C B(D° -> tf-7r+) + K 2 i3(D° -> A+tt-) - 2C 2 4gSi3(D+ -> A+tt ) 

cos 5^= , (26) 

2 tan 2 9 c kJb(D° -> K-rr+)B(D° -» A+tt-) 



where r(D) is the life time of D. Obviously, in the SU(3) symmetry limit, « = sc = 1, £ = 0, 
eq. (HH) will be A(D° -> K + tt~) = - tan 2 9 C A(D° -> i^"7r + ), and 5 x?r vanishes. Note that 
the relation (l24l) does not depend on the color-suppressed amplitudes C and C because these 
amplitudes have nothing to do with the above three decay modes. 

In order to go further into the analysis, we need to have information about the form 
factors F®^ p (q 2 ). For their q 2 dependence, we adopt the Bauer-Stech-Wirbel model [21], in 
which the form factors are assumed to behave as a monopole, 

Fr P (<f) = (27) 
1 — q /m* 

where m* is the pole mass with m* = 2.47 GeV for P = it and m* = 2.60 GeV for P = K. 
F®^ p (0) can be obtained via F^^ p (0) = F+~* p (0), since the latter can be measured in 
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-0.2 


-0.3 


-0.4 


-0.5 


-0.6 


COS 5 Kir 


0.979±0.018 


0.975±0.021 


0.969±0.025 


0.960±0.032 


0.946±0.043 


Kit 


11.7° ±4.9° 


12.9° ±5.3° 


14.3° ±5.8° 


16.2° ±6.6° 


18.9° ± 7.7° 


COS 8 Kir 


0.983±0.015 


0.980±0.017 


0.976±0.021 


0.970±0.026 


0.960±0.034 


&Kir 


10.5° ±4.6° 


11.4° ±4.9° 


12.6° ±5.5° 


14.1° ±6.1° 


16.2° ±6.9° 



Table 2: The relative strong phase between D° — > K + ir~ and D° —> K~tt + predicted by eq. 
(KoT) for different values of C 2 /Ci. B(D + — > K + n°) by BaBar [T5] is used for the results in 
the second line; B{D + — > K + tt°) by CLEO [H] is used for the results in the third line. The 
sign of 5 Kir could also be minus. 

semi-leptonic D° — > 7r _ £ + z/ and D° — > K~£ + v decays. The new experimental values from 
the Belle Collaboration [25] give 

F£ ~* k (0) = 0.695 ± 0.023, F%^*(0) = 0.624 ± 0.036, 

F + ~**(0)/F? ^ K (0) = 0.898 ± 0.045, (28) 

which are consistent with very recent results from lattice calculation |26j and from the QCD 
sum rules calculation [27] . In practice, only the ratio of these two form factors in ( 1281) is 
needed for our numerical calculations. By applying it to eq. (I25I) . we get 

x = 1.002 ±0.050, (29) 

which is very close to its value in the SU(3) symmetry limit [the error in (1291) is due to the 
uncertainty of F®^ n (0) / F®~* K (0) only]. Although this may be just a numerical coincidence, 
it seems that SU(3) symmetry breaking effects are dominated by decay constants (f n and 
fx). 

The numerical predictions of Sk-k for different values of GijC\ are displayed in Table 
2. As mentioned above, C2/C1 is regarded as a varying parameter. B(D° — > K~tt + ) and 
B(D° — > K + tt~) are taken from [7J. Since Particle data group has not given the average for 
B(D + — > K + n°) yet, both of the measurements listed in eq. (jl6"j) have been used, and the 
results are shown in the second and third lines of Table 2, respectively. The error is due to the 
uncertainty of x in eq. ( |29l) and the uncertainties of experimental branching ratios mentioned 
above, in which the contribution from B(D° — > K~n + ) and B(D° — > K + tt~) is actually very 
small and can be neglected. From Table 1, for the relevant scale of D decays, i.e. \x in the 
range of 1.0 ~1.5 GeV, CijG\ is about —0.5 ~ —0.3. Therefore, a not large but nonzero 
5 Kir, whose magnitude is 10° or above, i.e. sin<5 ~ ±0.2, might be expected from the present 
analysis. The authors of Ref . [1] , by assuming the existence of nearby resonances for the D 
meson, have obtained sin Sk-k — ±0.31, namely, cos 5ktt = 0.951 and Sk-k is about ±18°. The 
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other existing hadronic models which incorporate SU(3) symmetry breaking effects seem to 
prefer a small value of this phase, sin S < 0.2, with most models giving sin 5 < 0.1 [28j [29] 
(see Table I in Ref. [2H] for details). Unfortunately, the current measurement of 5k-k is very 
rough [13], as shown in eq. Employing the asymmetry R(D) measurements with some 
theoretical assumptions, another experimental result 5k-k ~ (3 ± 6 ± 7)° with relative small 
uncertainty is induced in Ref. |11J . Both of them are still consistent with zero. 

Finally, we estimate D — > Ktt decay amplitudes from the currently available data. The 
three independent complex phenomenological parameters are chosen as T, C, and £, not as 
a^ ff , a| ff , and £, because we will only use the ratio of the form factors [in eq fl28l) ] instead of 
their absolute values in the analysis. Without loss of generality, T is set to be real. 8q (Se) 
is the relative strong phase of C {£) to T. Here we take B(D° -> K'ir + ), B(D° -> K°7r°), 
B(D + — > K°7T + ), and B(D° — > K + n~) given by Particle data group jT], together with 
B(D + -> K + n°) by BaBar Collaboration [15J to illustrate our numerical calculation. The 
results of T, C, £, and |a| / a i ff I are summarized in Table 3, and other amplitudes T 7 , C, 
£', and A' can be easily derived using eqs. ( fT7l) . ( fT9l) and ( |2TI) . Several observations and 
remarks are given as follows. 



• The color-suppressed amplitude has a phase ~ 160° relative to the color-allowed ampli- 
tude T, and \C\ is effectively enhanced. This means that there could exist the strongly 
destructive interference between T and C. We get af^/a\ s ~ 0.56e ±ll60 °, which is 
insensitive to the value of C2/C1. a| ff /a^ ff = 0.62e~ il52 ° is obtained in [22]. 

• The £ amplitude has a relative phase ~ 130° to T, ~ 70° to C. Its magnitude is 
relatively large, and \£\ > \C\. This is contrary to the results in Ref. [22l [3j. As 
pointed out in [18], in general, the weak annihilation parameter Xa in eq. (120]) should 
be of order ln(m£)/A) and A is a soft scale. By taking a s ~ 0.5, C\ ~ 1.2, and 
£ = 0.325e ±il31 ° for C 2 jC x = -0.3, we can roughly estimate X A = 4.09e ±il22 ° or 
3.66e ±m °, which indicates \X A \ ~ 21n(m D /A) with A ~ 0.3 GeV. \X A \ = 3.84 has 
been obtained in [2"U] . 

• Some of our results are not in agreement with the ones in Refs. [221 E] , since we do 
not work in the SU(3) symmetry limit, and we mainly concentrate on D — > Krr decay 
modes in this paper. 



We return to discuss the asymmetry R(D + ). As pointed out before, the charged case is 
not as simple as the neutral case. Because of 

A \ D * - K Y\ = - tan* e c %±4 = - tan* C C + C ^ £ , (30) 
A{D+ -> ir°7r+) c C + T C + T ' v; 

one cannot simplify it as a similar analytic relation (T2"2"|) for neutral modes under C = C and 
£ = even if including the additional constraint f|2T|) already. However, using the values 
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Ca/d T [GeV 3 ] C [GeV 3 ] £ [GeV 3 ] |af /of | 

-0.3 0.417 0.289 e ±il60 ° 0.325 e Til31 ° 0.568 

0.445 0.306 e ±il63 ° 0.368 e Tll35 ° 0.563 

0.334 e ±il67 ° 0.425 e^ 139 ° 



-0.4 
-0.5 



0.485 



0.564 



Table 3: Numerical results of quark-diagram amplitudes T, C, S, and of | a| ff / ff | estimated 
by using the present data with different C 2 /Ci. Only the central values of the magnitude 
and the phase are quoted. 



of amplitudes listed in Table 3, numerically, we will get 

' -tan 2 c 1.538e ±il06 °, C 2 /Ci 



-tan 2 # c 1.532e ±il05 °, C 2 /d 



-0.3, 
-0.4, 



(31) 



which lead to 







tan 2 6 C 


1.521e ±il03 ° 


, C 2 /Ci = -0.5 






' 0.044, 


c 2 /c x = - 


-0.3, 


R(D+) 




0.040, 




•0.4, 






0.035, 


cyd = - 


-0.5. 



(32) 



The present observed value by CLEO Collaboration [II] is R(D+) = 0.030 ± 0.023 ± 0.025. 
Also, the suppression of R(D + ) comparing with R(D°) can be understood. From the defini- 
tion of R{D) in eq. (jlj), one will find it is proportional to 2 tan 2 9c cos 5, and S is the relative 
strong phase between the corresponding DCS amplitude and the CF amplitude. Now it is 
found that, 5 vanishes in the D shown in eq. fT22|) ; while 5 is about 100° in the D + 

case [see eq. (|3~T|) ]. Therefore R(D + ) is suppressed by small cos S. 

Furthermore, we discuss the possible generalization to the analysis of SCS D — > tttt, KK 
decays. Consider the ratio 

V rJ , 2 T(D+ 



Ri 



V, 



cd 



T(D- 



(33) 



The recent measurement gives Ri = 1.54 ± 0.27 [13 Ej and Ri should be unity in the SU(3) 
symmetry limit. Note that these two modes have only T and C amplitudes. From our 
analysis, one can get 



Ri 



1.073 



F°- 



m 



rn 



K ) 



X 



1 + af/af 



l/x+(f K /f n )af/af 



(34) 
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where the factor 1.073 is from the phase-space differences for the tttt and Kir final states. 
Taking a^/af 1 ~ 0.56e ±il60 ° from Table 3, we obtain Ri ~ 1.44, in accord with the recent 
measurement. Likewise, the ratio 

T(D° -> K+K~) , x 

can also be estimated using Table 3. Unfortunately, this result is far from the experimental 
value T(D° -> K+K-)/T(D° -> tt+tt-) = 2.82 ± 0.10 [7], implying that SU(3) symmetry 
breaking is still not fully accounted for. Since now there exist weak-annihilation contri- 
butions, in deriving eq. fl35l) . we have assumed £r+k- = fx/f^ and £ n7r = fn/fx^- 
This is actually not true because, under this assumption, the amplitude for the pure weak- 
annihilation D° -> K°K° decay will vanish, whereas B(D° -> = (7.4 ± 1.4) X 10~ 4 
experimentally [7J. Therefore the failure of reproducing the experimental value in eq. ( 1351) 
may be unavoidable in the present framework. In the case of R\, the weak-annihilation 
contribution is fortunately absent. This implies that the above relations for £k+k- and £ 7T7T 
need some corrections, and the weak-annihilation amplitudes should be carefully investigated 
when one would like to generalize the present work to the case of the SCS D — ► nn, KK 
decays including the pure weak- annihilation mode D° — > K°K°. 



4 Summary 

We have presented a phenomenological analysis of D — ► Kn decays including both CF and 
DCS modes. In order to determine all decay amplitudes for these processes using the present 
data, a moderate SU(3) symmetry breaking formalism has been assumed. Our analysis 
indicates this assumption works well in D — > Kn decays. The color-suppressed amplitude 
is enhanced, and it has a phase ~ 160° relative to the color-allowed amplitude. A large 
weak annihilation amplitude is obtained. Both of the asymmetries R(D°) and R(D + ) have 
been predicted, which are in good agreement with the experimental data. Our analysis also 
shows that a not large but nonzero 8 Km which is about 10° or above, might be expected. 
This means that there is no good reason to take sin 5k-k = in the experimental analysis of 
D° — > K ±r n^ decays. A precise measurement of 8kk will be welcome both theoretically and 
experimentally. 

We would like to point out that, the relation (1211) between £' and A' amplitudes, which 
is important to enable us to calculate 5k-k and estimate the D — > Ktt amplitudes in this 
analysis, is a model dependent assumption. Further tests for this relation will be very useful. 
But such tests cannot be performed at present because of the lack of suitable data. On the 
other hand, a similar relation between lU-exchange and PU-annihilation amplitudes could 
occur in SCS D° —>■ K + K~ and D + — > K + K° decays, which might provide an interesting 
test. As mentioned above, the present analysis however cannot be generalized to SCS decays 
straightforwardly when these transitions receive contributions from the weak-annihilation 
amplitudes. Therefore it would be useful to extend the present framework assuming (T2T1) 
to include also SCS decays D — > tttt, KK. Since here our main analysis concerns D — * Ktc 
decays, a further discussion of these issues will be left for future work. 
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Very recently, the similar study for R(D°) has been obtained in Ref. [23]. Since in 
our framework, we can employ the currently available data to determine all quark-diagram 
amplitudes T(T'), C(C), £(£'), and A(A') including their relative strong phases, R(D + ) 
and Sk-k are also estimated. 
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